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Recent Results on Flavor Physics from BaBAR 

J. Benitez on behalf of the BaBAR Collaboration 
SLAC, Stanford, CA 94025, USA 

We report an update to our previous measurement of the CKM element \Vub\ using exclusive 
B — >■ ttIu decays. In the charm sector we have performed a measurement of fo^ using Df — > t^Vt 
decays, we have measured the mixing parameter ycp using the lifetime ratio ^^^'^^ in -D" decays, 
and we have also searched for CP violation using T-odd correlations in D" decays to K'tt^tt' . 
Finally, in the tau sector we have performed a search for the lepton flavor violating decays — >■ e*7 

, I. INTRODUCTION 

o 

(N 

^ ■ The physics reach of the EaBAR experiment encompasses a large part of the flavor sector of particle physics. 
^ The BaBAR detector at the SLAC PEP-II asymmetric-energy e+e~ collider collected approximately 500 fb~^ at 
center of mass (CM) energies near 10.58 GeV between 1999 and 2008. At this energy the production cross section for 
production of fermion-anti-fermion pairs cc, bb and t+t~ is roughly the same and yields an event sample of about 600 
^—1' Million produced events of each type. These samples enable precise measurements of Cabibbo-Kobayashi-Maskawa 
^ : (CKM) parameters, studies of charm meson properties, and searches for rare r decays among other topics. This 
article presents a summary of recent measurements in these areas from BaBAR. 
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II. DETERMINATION OF IKJ FROM B 7r£iy DECAYS 



' The elements of the CKM quark-mixing matrix are fundamental parameters of the Standard Model (SM) of 
^ ^ ' electroweak interactions. With the increasingly precise measurements of decay-time-dependent CP asymmetries in 
_B-meson decays, in particular sin(2/3) [HQ, improved measurements of \Vub\ and |ycfc| will allow for more stringent 
] experimental tests of the SM mechanism for CP violation The best method to determine \Vub\ is to measure 
■ semileptonic decay rates for B — X^d^ {X^ refers to hadronic states without charm), which is proportional to |KihP- 
We have performed a study of four exclusive semileptonic decay modes, B'^ Tr~£~^v, ir'^t^h', 5° -> p~£'^i^, 

,—1 , and _B+ p^t'^v, and a determination of \Vub\- Here I refers to a charged lepton, either e+ or and v refers to the 
^ ■ associated neutrino. Exclusive decays offer good kinematic constraints, and thus effective background suppression. 
This analysis represents an update of an earlier measurement 4] that was based on a significantly smaller data 
set. For the current analysis, the signal yields and background suppression have been improved, and the systematic 
5^ , uncertainties have been reduced through the use of improved reconstruction and signal extraction methods. 

In the analysis of B^ — > tt^ i'^v we reconstruct the pion and lepton tracks, and determine the neutrino 4-momentum 
as the missing momentum in the event. A similar reconstruction is performed for the other channels. We then 
determine the signal yield as a function of by performing a two-dimensional fit in the variables toes and Ai?. 
AE is the difference between the reconstructed B energy and half the CM beam energy, and toes is the mass of 
the B candidate computed using the reconstructed 3-monientuni and half the beam energy. The projections of the 
fit for B^ -> n~l'^v are shown in Fig[TJ The fit is performed simultaneously to all channels while constraining 
B+ — >■ TT^^+i/ using isospin symmetry. The B^ — >■ p~t'^v and B'^ — >■ p^t^v modes help constrain the background. 
From the extracted signal yields the decay rate can be determined using the known total number B events produced 
and correcting for the reconstruction efficiency. 

The decay rate for B'^ tt^£^i' depends on the momentum, q^, carried by the and takes the following formula 
in the SM, 

dq^d cos 9wi ^2^^'^'' ^' ' 
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FIG. 1: Left plots show the rriES distributions in each bin for B*^ — > tv' i^v. The points show the data and the histograms 
show the fit and the different background contributions. The right-most plot shows the simultaneous fit to data (points) and 
to the FNAL/MILC lattice prediction (magenta, closed triangles). The LQCD results are rescaled according to the |Ki6| value 
obtained from the fit. 



where is the momentum of the pion in the rest frame of the B meson and Owi is the angle of the charged-lepton 
momentum in the W rest frame with respect to the direction of the W boost from the B rest frame. The g^-dependent 
form factor, f+{q^), is calculated from Lattice QCD. Using this formula we convert the form factors into decay rates 
and require that they match the measured values as shown in FigH] From the conversion factor \Vub\ is determined 
to be (2.95 ±0.31) x IQ-^. 

III. MEASUREMENT OF THE BRANCHING FRACTION FOR £)+ t+u^ AND EXTRACTION OF 

THE DECAY CONSTANT Fds 

The purely leptonic decays of the 13^ meson provide a clean probe of the pseudoscalar meson decay constant , 
which describes the amplitude for the c and s quarks to have zero spatial separation within the meson. In the SM 
these decays occur through a virtual boson which decays to a lepton pair, ignoring radiative processes, the total 
width is 

where and mt are the and lepton masses, respectively, Gf is the Fermi coupling constant, \Vcs\ is the 

magnitude of the CKM matrix element that characterizes the coupling of the weak charged current to the c and 
s quarks. In the context of the SM, predictions for meson decay constants can be obtained from QCD lattice 



calculations [6l-l9l. Illl|. The most precise theoretical prediction for fo^ is (241±3) MeV [8|. This value is in slight 



disagreement with the current measurement of Jd^ [l^ . It is important to validate the lattice QCD predictions 
through measurements of Jd^ ^ these computational methods are also used in other areas such as B meson decays. 
In addition, it is possible that physics beyond the SM can induce a difference between the theoretical prediction and 
the measured value. 

The relatively large branching fraction for the t"*" decay mode motivates the use of the decay sequence — >■ t'^Vt, 
T+ — e'^VeVr in this analysis. We use the well known branching fraction B{Df — ^ K'^K^) for normalization. This 
analysis uses an integrated luminosity of 427 fb^^ corresponding to the production of approximately 554 million cc 
events. 

Signal events are reconstructed in the production processes e+e" — > cc — > D*~^ -Dtag K^'~ X, with the subsequent 
decay — >■ -D+7. Here, I?tag is a fully reconstructed hadronic D meson decay, required to suppress the large 
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FIG. 2: Left two plots show the recoil mass for Df t'^Vt (r^ e^UeVr) with -Bextra ~ and -EcxtraO. -Eoxtra is the remaining 
energy in the calorimeter after the full event reconstruction. The two right plots show the corresponding distributions for 
Dt K^K+ events. The solid curve shows the total fit while the dashed-curve shows the signal component. 



background from non-charm continuum qq pair production; X represents a set of any number of pions (tt'^ and tt*) 
produced in the cc fragmentation process, and K^'^ represents a single K'^ or from cc fragmentation required 
to balance strangeness in the event. In addition we require a reconstructed vifhich tags the decay r+ — > e^Vet^r- 
In a similar reconstruction we select — > K'^K^ events. 

We extract the signal yields using the candidate mass determined from the 4-momentum recoiling against the 
D'TAcK'^' X-f system. The fit results are shown in Fig. [5J We compute the branching fraction using the formula, 

BjPt ^ r+j^r) _ BjKl ^ TT+n-) (NsY-^ e^sK+ 

B{Dt ^ K^,K+) B{t+ ^ e+y,Vr) (iVs)^g^+ e^'^^ ' ^ ' 

where Ns and e refer to the number of signal events and total efficiency for the tv and the normalizing decay 
modes. The values of the — >■ 7r"'"7r~ and r"*" e^vjjr branching fractions are obtained from [l^]. We find 
B(D^ T+Vr) ^ (4.5±0.5±0.4±0.3)% and use Eq.[2]to compute fo, =(233±13±10±7) MeV Here the errors 
are statistical, systematic, and due to PDG parameter values. 



IV. MEASUREMENT OF D^-D^ MIXING USING THE RATIO OF LIFETIMES 
FOR THE DECAYS K-tt+ AND K+K' 



Mixing in the charm sector has only recently been observed at B-factories. One manifestation of D'^-D^ mixing 
is differing decay time distributions for decays to different CP eigenstates 15|. We present here a measurement 
of this lifetime difference using a sample of and decays in which the initial flavor of the decaying meson is 
unknown. 

Assuming CP conservation in mixing, the two neutral D mass eigenstates and ID2) can be represented as 

I A) = pm+qm 

\D2) = p\D^)^q\D^) , ^ ' 

where \p\ + \q\ = 1. The rate of D^-D^ mixing can be characterized by the parameters x = Am/T and y = Ar/2r, 
where Am = mi — m2 and AT = Fi — r2 are respectively the differences between the mass and width eigenvalues 
of the states in Eq. and F = (Fi + F2)/2 is the average width. If either cc or y is non-zero, mixing will occur, 
altering the decay time distribution of and D° mesons decaying into final states of specific CP 13| . In the limit 
of small mixing, and no CP violation in mixing or in the interference between mixing and decay, the mean lifetimes 
of decays to a CP eigenstate of samples of D'^ {'''hh) Si.nd {Thh) mesons, and the mean lifetime of decays to a 
state of indefinite CP {tktt), can be combined to form the quantity 

ycp = 7 r - 1, (5) 
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FIG. 3: The left plots show the D candidate invariant mass distribution for K~n and K K~ , data are shown by points, 
total fit as by a curve and background contribution as solid shade. The right plots show the decay time distribution for data 
(points), total lifetime fit (curve), combinatorial background (gray), and charm background (black) contributions overlaid. 



where {rhh) = (^"/^ + '''/fh )/2- An analogous expression {jk-k) holds for the K^-k^ final state. If ycp is zero there is 
no D'^-D'^ mixing attributable to a width difference, although mixing caused by a mass difference may be present. 
In the limit of no direct CP violation, ycp = y- 

In this analysis we reconstruct mesons in reactions of the kind e+e^ — >■ cc — )■ D'^X where X is any additional 
system and decays to either K~K~^ or K~Tr'^. The decay time of the D'^ candidates is calculated using the 
measured displacement of the decay vertex with respect to the e~^e~ interaction region. The _D" candidate mass 
distributions and decay time distributions are shown in Fig. |3l To determine the Z?*^ lifetime in each channel we 
fit the decay time distributions using an exponential distribution convolved with the resolution function determined 
from simulated signal events. We find a value of j/cp (untagged) = [1.12 ± 0.26(stat) ± 0.22(syst)]%, which excludes 
the no-mixing hypothesis at 3.3a [14 1. 



V. SEARCH FOR CP VIOLATION USING T-ODD CORRELATIONS IN 

D° ^ K+K-7V+7V- DECAYS 

Physics beyond the SM, often referred to as New Physics (NP), can manifest itself through the production of 
new particles, probably at high mass, or through rare processes not consistent with SM origins. SM predictions 
for CP asymmetries in charm meson decays are generally of O(10~^), i.e. at least one order of magnitude lower 



than current experimental limits 27[. Thus, the observation of CP violation with current sensitivities signal NP. 
We report the results of a search for CP violation in the decay process — K~^K~tt^tt~ using a kinematic triple 
product correlation of the form Ct — Pi ■ (P2 x Pa), where each pi is a momentum vector of one of the particles in 
the decay. The product is odd under time-reversal (T) and, assuming the CPT theorem, T-violation is a signal for 
CP- violation. Strong interaction dynamics can produce a non-zero value of the At asymmetries, 

^ TjCr > 0) ~ nCr < 0) 
^ ~ T{Ct > 0) + F(C7t < 0) ' ^ ' 

where F is the decay rate for the process, even if the weak phases are zero. After defining a similar formula for the 
CP-conjugate decay process we can construct At = ^{^t — At), a non-zero value of At would signal CP-violation 
0. 

Following the suggestion by I.I. Bigi [18] to study CP violation using this technique, the FOCUS collaboration 
made the first measurements using approximately 800 events and reported At{D° K^K^tt^tt^) = 0.010±0.057± 
0.037 [l^. We perform a similar study using approximately 1.5 x 10^ signal events. 

Reactions of the kind e+e~ -> X D*+; £>*+ 7r+L)°; £)° K+K^n+n' , where X indicates any additional 
(unreconstructed) system, have been selected. We require the D'^ to have a CM momentum greater than 2.5GeV/c. 
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FIG. 4: Fit projections onto in{K'^ tv^tt^) for the four different Ct subsamples after a Am signal selection. The shaded 
areas indicate the total backgrounds. The normalized fit residuals, represented by the pulls are also shown under each 
distribution. 



According to the Z3*+ tag and the Ct variable, we divide the total data sample into four subsamples. The D° yields 
are determined using a binned, extended maximum-likelihood fit to the 2-D {m{K'^ K^tt'^tt^), Am) distribution 
obtained with the two observables m{K'^ K~Tr'^7T~) and Am = m{K~^ K~ tt"*" ■n~ 7r+) - m{K^ K~'K'^Tr~). The 
functional forms of the probability density functions (PDFs) for the signal and background components are based 
on studies of MC samples. We make use of combinations of Gaussian and Johnson SU [i^l line shapes for peaking 
distributions, and we use polynomials and threshold functions for the non-peaking backgrounds. Fig. |4] shows the 
K'^ K~ TF"*" 7r~ mass distributions for the four different Ct subsamples. The asymmetry determined using the signal 
yields is found to be consistent with zero: At = (1.0 ± 5.1g^a,t ± 4.4gyg^) x 10"'^ with a sensitivity of ~ 0.6% |16|. 



VI. SEARCHES FOR LEPTON FLAVOR VIOLATION IN THE DECAYS 

e=^7 AND ^i^j 

Despite the existence of neutrino oscillations [22I1, decays of — ?> £^7 (where £ — e,^) are predicted to have 
unobservably low rates [231] in the SM. Thus, observation of charged lepton flavor violation would be an unambiguous 
signature of new physics. Presently, the most stringent limits are 6(r± ^ e±7) < 1.1 X 10-^ 24,] and S(r± 
/i±7) < 4.5 X 10-8 Q at 90% C.L., using 232.2 fb"^ and 535 fb"^ of e+e" annihilation data collected near the 
T{AS) resonance by the BaBAR and Belle experiments, respectively. This analysis utilizes the entire BaBAR dataset 
corresponding to a luminosity of 425.5 fh'^ , 28.0 fb"^ and 13.6 fb"^ recorded at the r(45), r{3S) and T{2S) 
resonances, and 44.4 fb^^, 2.6 fb^^ and 1.4 fb^^ recorded at 40MeV, 30MeV and 30MeV below the resonances, 
respectively. 

The signal is characterized by a pair with an invariant mass and total energy in the CM frame (Ef^) close to 
rrir = 1.777 GeV/c^ [13] and •\/s/2, respectively. Candidate events must also contain another t decay product (one or 
three tracks). The signal-side hemisphere must contain one photon with CM energy greater than 1 GeV and 

one track within the calorimeter acceptance with momentum in the CM frame less than 0.77^/s/2. This track must 
be identified as an electron or a muon for the — > e^7 or — /z^7 search. 

Signal decays are searched for using two kinematic variables: the energy difference AE = Ef^ — y/sjl and the 
beam-energy constrained t mass (toec), obtained from a kinematic fit after requiring the CM r energy to be ^fsjl. 
For signal events, the rriEc and AE distributions are centered at to^ and small negative values, respectively, where 
the shifts from zero for the latter are due to radiation and photon energy reconstruction effects. The toec vs. AE 
distributions are modeled by 2-dimensional probability density functions (PDFs) summed over all background event 
types. We observe and 2 events for the — )■ 6^7 and — )■ /i^7 searches inside the 2a signal ellipse as shown in 
Fig. [S] As there is no evidence for a signal, we set a frequentist upper limit calculated using B'jji^ — N^^/{Nre) to 
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FIG. 5: Distributions of 6*7 (left) and (rigiit) candidate decays in the rriEC vs. AE plane. Data are shown 

as dots and contours containing 90% (50%) of signal MC events are shown as light- (dark-) shaded regions. The 2a ellipse is 
shown also. 



be S(t± -> e±7) < 3.3 x IQ-^ and S(r± ^±7) < 4.4 x IQ-^ at 90% C.L. fill, where e is the signal efficiency 
inside the 2a signal ellipse and Nfjj^ is the 90% C.L. upper limit on the number of signal events, estimated using the 
POLE program 0. 



VII. CONCLUSIONS 



In conclusion, the BaBAR collaboration continues to exploit its rich data-set to study fundamental aspects of flavor 
physics. In this article we report an update to our previous measurement of the CKM element \Vub\ using exclusive 
B Titv decays. In the charm sector we have extracted a value of fo^ from — >■ r+^'^- decays, we have measured 
the mixing parameter ycp using the lifetime ratio -^^^^ in -D" decays, and we have also searched for CP violation 
using T-odd correlations in 4-body decays to tt^-k^ . Finally, in the tau sector we have placed upper limits 

on the rates of lepton flavor violating decays — > 6^7 and — > /i^7. 
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